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(Ga,Mn)As is at the forefront of research exploring the synergy of magnetism 
with the physics and technology of semiconductors, and has led to discover- 
ies of new spin-dependent phenomena and functionalities applicable to a wide 
range of material systems. Its recognition and utility as an ideal model material 
for spintronics research has been undermined by the large scatter in reported 
semiconducting doping trends and micromagnetic parameters. In this paper 
we establish these basic material characteristics by individually optimizing the 
highly non-equilibrium synthesis for each Mn-doping level and by simultane- 
ously determining all micromagnetic parameters from one set of magneto-optical 
pump-and-probe measurements. Our (Ga,Mn)As thin-film epilayers, spannig 
the wide range of accessible dopings, have sharp thermodynamic Curie point 
singularities typical of uniform magnetic systems. The materials show system- 
atic trends of increasing magnetization, carrier density, and Curie tempera- 
ture (reaching 188 K) with increasing doping, and monotonous doping depen- 
dence of the Gilbert damping constant of ~ 0.1 — 0.01 and the spin stiff'ness of 
~ 2 — 3 meVnm^. These results render (Ga,Mn)As well controlled degenerate 
semiconductor with basic magnetic characteristics comparable to common band 
fer r omagnets . 

Under equilibrium growth conditions the incorporation of magnetic Mn ions into III- 
As semiconductor crystals is limited to approximately 0.1%. To circumvent the solubility 
problem a non-equilibrium, low-temperature molecular-beam-epitaxy (LT-MBE) technique 
was employed which led to first successful growths of (In,Mn)As and (Ga,Mn)As ternary 
alloys with more than 1% Mn and to the discovery of ferromagnetism in these materials.'^ 
The compounds qualify as ferromagnetic semiconductors to the extent that their magnetic 
properties can be altered by the usual semiconductor electronics engineering variables, such 
as doping, electric fields,'^'^ or light .'^^HSI gy exploiting the large spin polarization of car- 
riers and low saturation moment in (Ga,Mn)As and building on the well established het- 
erostructure growth and microfabrication techniques in III-V semiconductors, (Ga,Mn)As 
has been extensively used for spintronics research of direct and inverse magneto-transport 
phenomena.'^^'^ Besides the more conventional spintronic effects based on Mott's two-spin- 
channel model of conduction in ferromagnets, (Ga,Mn)As has become particularly fruitful 
for exploring the second, more physically intriguing spintronics paradigm based on Dirac's 
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spin-orbit coupling .I^SESHIS] 

The apparent potential of (Ga,Mn)As to become the test-bed model material for many 
lines of spintronics research has been hindered by the large scatter in reported semiconduct- 
ing doping trends and micromagnetic parameters. Our strategy to tackle this problem begins 
from the synthesis of a set of (Ga,Mn)As materials spanning a wide range of Mn dopings. Be- 
cause of the highly non-equilibrium nature of the heavily-doped ferromagnetic (Ga,Mn)As, 
the growth and post-growth annealing procedures have to be individually optimized for each 
Mn-doping level in order to obtain films which are as close as possible to idealized uniform 
(Ga,Mn)As mixed crystals with the minimal density of compensating and other uninten- 
tional defects. An extensive set of characterization measurements has to accompany the 
synthesis to guarantee that the materials show systematic doping trends; monitoring the 
thermodynamic Curie point singularities is essential for assuring the uniformity and high 
magnetic quality of the materials.'^^SZHlS] wj^en omitting the above procedures,'^ extrinsic 
impurities and sample inhomogeneities can yield non-systematic doping trends and conceal 
the intrinsic magnetic properties of (Ga,Mn)As. 

The focus of the work presented in this paper is on the systematic study of the Gilbert 
damping and spin stiffness constants which, together with magnetic anisotropy fields, rep- 
resent the basic micromagnetic parameters of a ferromagnet. A more than an order of mag- 
nitude experimental scatter and a lack of any clear trend as a function of Mn-doping can be 
found in the literature for the Gilbert damping and spin-stiffness constants.'^^'^ (See Supple- 
mentary information for a detailed discussion of previous experimental works.) This reflects 
partly the issues related to the control and reproducibility of the synthesis of (Ga,Mn)As 
and partly the difficulty with applying common magnetic characterization techniques, such 
as neutron scattering, to the thin-film dilute-moment (Ga,Mn)As samples. Hand-in-hand 
with the optimization of the material synthesis we have developed experimental capabilities 
based on the magneto-optical (MO) pump-and-probe method which allow us to simultane- 
ously determine the magnetic anisotropy, Gilbert damping, and spin stiffness constants from 
one consistent set of measured data. Our results are summarized in Fig. 1. The Curie point 
singularity in the temperature derivative of the resistivity dp/dT measured throughout the 
series of optimized ferromagnetic (Ga,Mn)As samples with metallic conduction is shown in 
Fig. la. The data span the nominal doping range from x ?«1.5 to 13% and corresponding 
Curie temperatures from T^. = 29 to 188 K, and illustrate the high quality of all the epilayers 
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FIG. 1: Micromagnetic parameters of optimized epilayers of ferromagnetic (Ga,Mn) As. 

a, Examples of sharp Curie point singularities in the temperature derivative of the resistivity in the 
series of optimized ferromagnetic (Ga,Mn)As epilayers with metallic conduction; Tc monotonously 
increases with increasing nominal Mn doping between 1.5 and 13%. Inset shows dp/dT normalized 
to its peak value with the temperature axis normalized to Tc. b, Examples of oscillatory parts of 
MO signals measured in 18 nm thick (Ga,Mn)As epilayers with the depicted nominal Mn doping 
for external magnetic field poH^xt = 400 mT applied along the [010] crystallographic direction; 
the curves are normalized and vertically off-set for clarity, c, Dependence of anisotropy constants 
on nominal Mn doping, d, Dependence of the Gilbert damping constant a and the spin stiffness 
constant D on nominal Mn doping. 



within the series. Examples of the measured magnetization precession signals by the MO 
pump-and-probe method are shown in Fig. lb. From these time-dependent magnetization 
measurements we obtained the magnetic anisotropy constants i^j, Gilbert damping constant 
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a, and spin stiffness constant D which are summarized in Figs.lc,d. We now proceed to the 
detail discussion of our experimental techniques and the discussion of the measured results 
in the context of physics of degenerate semiconductors and band ferromagnets. 

Optimization of the (Ga,Mn)As synthesis. Our (Ga,Mn)As layers were grown at 
the growth rate of approximately 0.2 monolayers/second. The Mn flux, and hence the nom- 
inal Mn doping x, was determined by measuring the ratio of the beam equivalent pressures 
(BEP) of Mn and Ga sources before each growth. The Mn content was cross-checked by 
secondary ion mass spectroscopy (SIMS) and by comparing the growth rates of GaAs and 
(Ga,Mn)As measured by the oscillations of the reflection high-energy electron diffraction 
(RHEED). 

There are two critical growth parameters of (Ga,Mn)As: the substrate temperature, 
and the As-to-(Ga-l-Mn) flux ratio. At the typical temperatures of ~ 200°C neither an 
optical pyrometer nor a radiatively coupled temperature sensor are applicable. Instead, 
we used the GaAs band-edge spectrometer to measure the substrate temperature and the 
predictive substrate heater control to stabilize the temperature during the growth. For a 
given As:(Ga-|-Mn) ratio the substrate temperature fully determines the growth regime: the 
growth proceeds two-dimensionally at low temperatures, and turns irreversibly into the 3D 
growth mode when a critical temperature is exceeded. The scatter of the critical substrate 
temperature for given x and As:(Ga-|-Mn) ratio is remarkably small, typically less than 
2°C. In excess As flux the 2D/3D transition occurs at higher temperature. The highest 
quality samples are grown in a narrow window of the 1:1 stoichiometric As:(Ga-|-Mn) ratio 
and at the substrate temperature approaching as close as possible from below the 2D/3D 
critical temperature for given x. The As:(Ga+Mn) ratio was adjusted by the As-cell valve, 
and calibrated using the As-controlled RHEED oscillations. In insets of Fig. 2a we show 
examples of RHEED patterns for the x — 7% nominally doped (Ga,Mn)As material grown 
at stoichiometric 1:1 ratio of As:(Ga-|-Mn) for substrate temperature of 225 K which is 
above the 3D/2D boundary and 210 K which is below the boundary. The optimal growth 
temperature for this doping is 215 K. In the main panel of Fig. 2a we plot the optimal 
growth temperature as a function of nominal Mn doping, showing the rapidly decreasing 
growth temperature trend. 

The next important factor determining the quality of the resulting (Ga,Mn)As materials 
are post-growth annealing conditions. In Fig. 2b we show the dependence of the Curie tem- 
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FIG. 2: Optimization of the (Ga,Mn)As synthesis, a, Optimal growth temperature as a 
function of the nominal Mn doping. Insets show examples of RHEED images of the 2D growth 
at 210°C (lower inset) and 3D growth at 225°C (upper inset) of the 7% Mn-doped (Ga,Mn)As. 
b, Dependence of the Curie temperature on the annealing time for three different annealing tem- 
peratures in a 15 nm thick (Ga,Mn)As epilayer with 13% nominal Mn doping, c, Dependence of 
the Curie temperature on the annealing time for the annealing temperature of 200° C in a 100 nm 
thick (Ga,Mn)As epilayer with 13% nominal Mn doping, and in the same epilayer thinned down to 
25 nm by wet etching, d— f, Temperature dependencies of resistivity p, d, temperature derivative 
of the resistivity dp/dT, e, and remnant magnetization M, f, in a 20 nm thick (Ga,Mn)As epilayer 
with 13% nominal Mn doping at successive annealing times at the optimal annealing temperature 
of 160° C for this doping. 

perature on the annealing time for three different annealing temperatures for the record 
Tc = 188 K sample with nominal 13% Mn doping and film thickness 15 nm. These curves 
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illustrate the common trend in annealing (at temperatures close to the growth temperature) 
suggesting the presence of competing mechanisms. One mechanism yields the increase of Tc 
and is ascribed in a number of reports to the removal of charge and moment compensating 
interstitial Mn impurities (see e.g. the detailed annealing study in Ref. [TOl) . The removal is 
slowed down by the growth of an oxide surface layer during annealin^^ and an additional 
mechanism can eventually yield reduction of Tc after sufficiently long annealing times, de- 
pending on the annealing temperature. The origin of this detrimental mechanism may be in 
Mn clustering or in the competition between the non-equilibrium (Ga,Mn)As phase and the 
equilibrium MnAs second phase. Because of the competing mechanisms, the absolutely high- 
est Curie temperature for the given nominal doping is achieved at intermediate annealing 
temperature and time, as illustrated in Fig. 2b. 

The remaining critical parameter of the synthesis is the epilayer thickness. For a given 
nominal doping, the highest attainable Tc is reached only in thin films, typically thinner 
than ~50 nm. In Fig. 2c we illustrate the importance of the film thickness for obtaining 
high quality (Ga,Mn)As materials. A 100 nm thick film is grown with nominal 13% doping 
and, unlike the thin record Tc film discussed above, here the maximum Tc achieved by 
annealing is only about 140 K. However, if the same film is thinned down (to e.g. 25 nm) 
by wet etching and annealed at the same conditions, the achieved Curie temperatures are 
significantly higher. 

An increase of Tc is not the only parameter followed to ascertain that a sample is of 
high quality. A key characterization tool are the thermodynamic Curie point singularities.'^ 
This is illustrated in Figs. 2d-f where we compare resistivity and magnetization measured 
at increasing time steps during the optimizing annealing procedure. The development of 
sharply vanishing magnetization M(T) at Tc and the onset of the singularity in dp/dT are 
well correlated with increasing Tc and conductivity within the annealing sequence. 

After finding the optimal growth and post-growth conditions for each individual nominal 
doping we obtained a series of samples spanning the wide range of Mn dopings. The samples 
can be divided into several groups: at nominal dopings below ~ 0.1% the (Ga,Mn)As ma- 
terials are paramagnetic, strongly insulating, showing signatures of the activated transport 
corresponding to valence band - impurity band transitions at intermediate temperatures, 
and valence band - conduction band transitions at high temperatures (see Fig. 3a) .'^^'^ For 
higher nominal dopings, 0.5 ^ x < 1.5%, no clear signatures of activation from the valence 
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FIG. 3: Doping trends in the series of optimized (Ga,Mn)As epilayers. a, Temperature 
dependence of the conductivity o"(T) of optimized (Ga,Mn)As epilayers with depicted nominal 
Mn doping. Daslied lines indicate the activated parts of cr{T) of the insulating paramagnetic 
(Ga,Mn)As with 0.05% Mn doping, corresponding to the Mn acceptor level and the band gap, 
respectively, b-e. Conductivity, b, hole density, c, saturation magnetization and corresponding 
Mn moment density, d, and Curie temperature, e, as a function of the nominal Mn doping in the 
series of optimized (Ga,Mn)As epilayers. 

band to the impurity band are seen in the dc transport, confirming that the bands start to 
overlap and mix, yet the materials remain insulating.'^^^ At x ~ 1.5%, the low-temperature 
conductivity of the film increases abruptly by several orders of magnitude (see Fig. 3b), and 
the system turns into a degenerate semiconductor.'^!'^ The onset of ferromagnetism occurs 
already on the insulating side of the transition at x ~ 1% and the Curie temperature then 
steadily increases with increasing nominal Mn doping up to ~ 13%. The hole concentration 
p can be measured by the slope of the Hall curve at high fields (see Supplementary infor- 
mation) with an error bar due to the multi-band nature estimated to ~ 20%.!^ Within this 
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uncertainty, the overall trend shows increasing p with increasing doping in the optimized 
materials, as shown in Fig. 3c. Similarly, the saturation moment and steadily increase 
with increasing nominal doping up to x ~ 13%, as shown in Figs. 3d,e. Assuming 
per Mn atoni^ we can estimate the density Nmu of uncompesated Mnca moments from the 
magnetization data (see left y-axis in Fig. 3d). An important conclusion can be drawn when 
comparing this estimate with the hole density estimated from the Hall resistance. Since 
there is no apparent deficit of p compared to Nmu, and since the interstitial Mn impurity 
compensates one local moment but two holes we conclude that interstitial Mn is completely 
(within the experimental scatter) removed in our optimally annealed epilayers. Hence, our 
series of optimized (Ga,Mn)As materials have reproducible characteristics, showing an over- 
all trend of increasing saturation moment with increasing x , increasing Tc (reaching 188 K), 
and increasing hole density. The materials have no measurable charge or moment compen- 
sation of the substitutional Mnca impurities and have a large degree of uniformity reflected 
by sharp Curie point singularities. 

Determination of the micromagnetic parameters. We now proceed to the de- 
termination of the magnetic anisotropy, Gilbert damping, and spin stiffness constants of 
our (Ga,Mn)As epilayers from the MO time-resolved measurements of the magnetization 
precession. In the MO pump-and-probe experiments, we used a femtosecond titan sapphire 
laser that was spectrally tuned to 1.64 eV, i.e., above the band gap of GaAs. The possibility 
to excite and detect precession of ferromagnetic Mn moments in (Ga,Mn)As by this method 
has been extensively discussed in previous MO studies.'^^^ All experiments presented be- 
low were preformed at temperature of approximately 15 K in reflection geometry. External 
magnetic fields up to 550 mT were applied in the [010] and [110] crystallographic directions. 
The intensity of the pump pulse was ~ 30 — 40 yuJcm"^, with the pump to probe intensity 
ratio ~ 20 : 1 — 10 : 1. The penetration depth of the laser beam (~ 600 nm) safely exceeds 
the thickness of the studied (Ga,Mn)As epilayers. 

The anisotropy constants, shown in Fig. Ic, where obtained combining three complemen- 
tary measurements. In the first experiment we measured the external magnetic field Hext 
dependence of the precession frequency / of the time resolved MO signal. In the studied 
(Ga,Mn)As/GaAs epilayers, the internal magnetic anisotropy fields are dominated by three 
components. The out-of-plane component Kout is a sum of the thin-film shape anisotropy 
and the magnetocrystalline anisotropy due to the compressive growth strain in (Ga,Mn)As. 
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The cubic magnetocrystalline anisotropy Kc reflects the zinc-blende crystal structure of the 
host semiconductor. The additional uniaxial anisotropy component along the in-plane di- 
agonal Ku is not associated with any measurable macroscopic strain in the epilayer and is 
likely of extrinsic origin. The precession frequency is given by, 

/ = ^^{H,,t cos((^ - ipn) - 2Kout + ^c(3 + cos4(^)/2 + 2Ku sm''{ip - 7r/4) + Aif„) 
X a/ {Hext cos(v? - (Ph) + 2Kc COS Aip - 2Ku sin 2(p + A//„) , (1) 

where g is the Lande g-factor of Mn moments, the Bohr magneton, ip and (pu are the 
in-plane magnetization and external magnetic field angles measured from the [100] crystal 
axis, and Ai7„ is the shift of the resonant field for the higher index n spin wave modes 
with respect to the n = uniform precession mode. In order to uniquely determine the 
anisotropy constants, the field-dependent precession frequency measurements were comple- 
mented by MO experiments with variable polarization angle of the probe beam. The latter 
measurements allow us to precisely determine the angle of the equilibrium easy axis of the 
magnetization (see Supplementary information) .1^^*^ Finally, we confirmed the consistency 
of the obtained anisotropy constants by performing static measurements of magnetization 
hysteresis loops by the superconducting quantum interference device (SQUID). Results for 
ferromagnetic materials from our series of optimized (Ga,Mn)As epilayers are summarized in 
Fig. Ic. Note that the values of Kout and for the given Mn-doping are well reproducible in 
materials whose synthesis yields the same optimized values of the basic structural, magnetic 
and transport properties. For the constant, variations in the width of the optimized 
thin (Ga,Mn)As films or of other otherwise insignificant changes of the growth or annealing 
conditions may yield sizable changes of K^. This confirms the presumed subtle extrinsic 
nature of this magnetic anisotropy component. 

The sign of Kout implies that all studied (Ga,Mn)As/GaAs materials are in-plane fer- 
romagnets. The competing magnitudes of Kc and Ku and the different doping trends of 
these two in-plane magnetic anisotropy constants (see Fig. Ic) are therefore crucial for the 
micromagnetics of the materials. The biaxial anisotropy K^. dominates at very low dopings 
and the easy axis aligns with the main crystal axis [100] or [010]. At intermediate dopings, 
the uniaxial anisotropy Ku is still weaker but comparable in magnitude to K^. In these 
samples the two equilibrium easy- axes are tilted towards the [110] direction and their angle 
is sensitive to small changes of external parameters such as temperature. This allows for 
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exciting the magnetization precession by laser pulses in the pump-and-probe MO experi- 
ments. At very high dopings, the uniaxial anisotropy dominates and the system has one 
strong easy-axis along the [110] in-plane diagonal. In the low-doped and high-doped samples 
with very stable easy-axes aligned with one of the main crystal directions the dynamical MO 
experiments become unfeasible. 

The Gilbert damping constant a, shown in Fig. Id, is obtained by fitting the measured 
dynamical MO signal to Landau-Lifshitz-Gilbert (LLG) equations using the experimentally 
obtained magnetic anisotropy constants. The high accuracy of the LLG fits is demonstrated 
in Figs. 4a,b on data measured in a x = 5.2% doped sample. The obtained dependence of a 
on the external magnetic field applied along the [010] and [110] directions is shown in Fig. 4c. 
At smaller fields, a is not constant and shows a strong anisotropy with respect to the field 
angle. When plotted as a function of frequency, however, the dependence on the field-angle 
disappears, as shown in Fig. 4d. Analogous results are obtained for the entire series of the 
optimized materials. We can therefore conclude that the apparent anisotropy of a can in 
our materials be ascribed fully to the field-angle dependence via the precession frequency. In 
all our studied materials, the frequency-independent Gilbert damping constant is isotropic 
and can be accurately determined from MO data with precession frequencies / ^ 15 GHz. 

We point out that in ferromagnetic resonance (FMR) experiments, the measurement fre- 
quency was limited to two values, / = 9 and 35 GHz which even in the optimized (Ga,Mn)As 
materials is not sufficient to reliably separate the intrinsic Gilbert damping constant from 
the inhomogeneous broadening of the FMR line-width. The dynamical MO measurements, 
on the other hand, span a large enough range of frequencies and allow us to extract a con- 
sistent set of frequency-independent values of a for our series of optimized ferromagnetic 
(Ga,Mn)As materials. We find a systematic doping trend across the series in which the 
Gilbert constant decreases from ~ 0.1 to 0.01 when the nominal Mn doping increases from 
~ 2% to 5% and then remains nearly constant (see Fig. Id). The magnitudes of a and the 
doping dependence are consistent with Gilbert damping constants in conventional transi- 
tion metal ferromagnets. In metals, a typically increases with increasing resistivity and is 
enhanced in alloys with enhanced spin-orbit coupling.'^SHnZl Similarly, in our measurements 
in (Ga,Mn)As, the increase of a correlates with a sizable increase of the resistivity in the 
lower Mn-doped samples. Also, the spin-orbit coupling effects tend to be stronger in the 
lower doped samples with lower filling of the valence bands and with the carriers closer to 
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FIG. 4: Determination of the Gilbert damping constant from MO experiments. a,b, 

Oscillatory part of the MO signal (points) measured in a 18 nm thick epilayer with 5.2% nominal Mn 
doping for external magnetic field ^QH^xt = 100 mT applied along the crystallographic directions 
[010] and [110]; lines are fits by the LLC equation, c, Dependence of the Gilbert damping on 
external magnetic field applied along the [010] and [110] crystallographic directions, d, Dependence 
of the Gilbert damping on the precession frequency. 

the metal-insulator transition.'^ Theory ascribing magnetization relaxation to the kinetic- 
exchange coupling of Mn moments with holes residing in the disordered, exchange-split, and 
spin-orbit-coupled valence band of (Ga,Mn)As yields a comparable range of values of a as 
observed in our measurement s.l^ 

Similar to the Gilbert constant, there has been a large scattei^ in previous reports of 
experimental values of the spin-stiffness in (Ga,Mn)As inferred from FMR,'^^!^!] magneto- 
optical studies,'^ and from complementary static magnetization and domain structure 
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measurements.'^^'^ We attribute the lack of a consistent picture obtained from these measure- 
ments to sample inhomogeneities and extrinsic defects in the studied (Ga,Mn)As epilayers 
with thicknesses typically exceeding 100 nm and to experimental data which allowed only 
an indirect extraction of the spin stiffness constant. The MO pump-and-probe technique 
utilized in our work allows in principle for the direct measurement of the spin stiffness, 
however, one has to find the rather delicate balance between thin enough epilayers to avoid 
sample inhomogeneity and thick enough films allowing to observe the higher-index Kittel 
spin-wave mode^ of a uniform thin-film ferromagnet. For these modes, the spin- stiffness 
parameter D is directly obtained from the measured resonant fields, 

r)2 77-2 

AH^ = Ho-H^ = D — , (2) 

where L is the thickness of the ferromagnetic film. The MO pump-and-probe technique 
has the key advantage here that, unlike FMR, it is not limited to odd index spin wave 
modes.'^ The ability to excite and detect the n = 0, 1, and 2 resonances is essential for the 
observation of the Kittel modes in our optimized (Ga,Mn)As epilayers whose thickness is 
limited to ~ 50 nm. 

In Fig. 5a we show an example of the time dependent MO signal measured in a 48 nm 
thick optimized epilayer with 7% nominal Mn doping. Three spin wave resonances (SWRs) 
are identified in the sample with frequencies /o, /i, and /2, as shown in Figs. 5b, c. The asso- 
ciation of these SWRs with the Kittel modes, described by Eq. (|2]), is based on experiments 
shown in Figs. 5b-e. In Fig. 5c we plot the dependence of the three detected precession 
frequencies on the external magnetic field applied along the [010] and [110] crystal axes. 
At saturation fields, which for the 7% Mn-doped sample are > 70 mT, the equilibrium 
magnetization vector is aligned with i^ext and Eq. ([T]) with = can be used to fit the 
data. We emphasize that all six displayed dependences fn{,Hext) for n = 0, 1, and 2, and 
ifH = 45° and 90° can be accurately fitted by one set of magnetic anisotropy constants. We 
can therefore use Eq. ^ to convert the measured frequency spacing of individual SWRs 
to AHn- In Fig. 5d we show that AHn in our optimized epilayers is proportional to as 
expected for the Kittel modes in homogeneous films. 

The magnetic homogeneity and the applicability of Eq. ^ in our epilayers is further 
confirmed by the following experiments: We prepared three samples by etching the original 
48 nm thick (Ga,Mn)As film down to the thicknesses of 39, 29 and 15 nm, respectively. 
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FIG. 5: Determination of the spin stiffness constant from MO experiments, a, Oscillatory 
part of the MO signal (points) measured in a 48 nm thick epilayer with 7% nominal Mn doping 
for external magnetic field /xo-f^ext = 20 mT applied along the [010] crystallographic direction; 
line is a fit by a sum of three damped harmonic functions, b, Fourier spectra of oscillatory MO 
signals (points) measured for fioHext = 20 mT applied along the [010] crystallographic direction 
in samples prepared by etching from the 48 nm thick epilayer. The curves are labeled by the film 
thicknesses, normalized, and vertically off-set for clarity; lines are fits by a sum of Lorentzian peaks, 
c, Dependence of the measured precession frequency (points) on the magnetic field for two different 
orientations of the field in the 48 nm thick epilayer; lines are fits by Eq. (1). d, Dependence of the 
measured mode spacing on square of the mode number in the 48 nm thick epilayer. e, Dependence 
of the spacing between the two lowest modes (AHi) on the film thickness. Lines in d and e are 
fits by Eq. (2) with spin stiffness D = 2 A3 meVnm^. 
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As seen in Fig. 5b, the frequency /o is independent of the film thickness which confirms 
that it corresponds to the uniform precession mode and that the film is homogeneous, i.e., 
the magnetic anisotropy constants do not vary across the width of the (Ga,Mn)As epilayer. 
The spacing Aifi shown in Fig. 5e scales as L^^ and the values of D extracted from the 
n-dependence of the resonant field spacings in the L = 48 nm epilayer (see Fig. 5d) and from 
the L-dependence of Ai/i (see Fig. 5e) give the same D = 2.43 ± 0.15 meVnm^. Identical 
value of the spin stiffness was also obtained from measurements in an epilayer grown with 
the same doping and thickness of 18 nm in which we detected the frequencies /o and /i 
and applied Eq. ([2]). These measurements confirm the reliability of extracted values of the 
spin stiffness. We note that the SWR frequencies are determined with high accuracy in 
our measurements and that the indicated error bars in Fig. Id reflect the uncertainty of 
the film thickness. As shown in Fig. Id, we observe a consistent, weakly increasing trend 
in D with increasing doping and values of D between 2 and 3 meVnm^ in the studied 
ferromagnetic samples with nominal doping 3.8-9%. (Note that apart from the difficulty 
of exciting magnetization precession in the very low and high-doped samples with stable 
easy-axes, the measurements of D were unfeasible on the lower doping side of the series 
because of the increasing damping and the corresponding inability to detect the higher SWR 
modes.) Similar to the Gilbert damping constant, our measured spin stiffness constant in the 
optimized (Ga,Mn)As epilayers is comparable to the spin stiffness in conventional transition 
metal ferromagnets.'^ 

We remark, that we tested the inapplicability of the SWR experiments for the direct 
determination of the spin stiffness in thick non-uniform materials. In the Supplementary 
information we show measurements in ~ 500 nm thick as-grown and annealed samples with 
7% nominal Mn-doping. The Curie temperatures of ~ 60 and 90 K can be inferred only 
approximately from smeared out singularities in dp/dT and M{T) and are significantly 
smaller than in the thin optimized epilayers with the same nominal doping. The films 
are therefore clearly inhomogeneous and contain compensating defects. Because of the large 
thickness of the epilayers we observe up to five SWR modes, however, consistent with the 
inhomogeneous structure of the films, the corresponding AiJ„ do not show the quadratic 
scaling with n of the Kittel modes of Eq. (|2|. 

In the experiments discussed above we have established the systematic semiconducting 
doping trends and basic magnetic characteristics of epilayers which have been optimized to 
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represent as close as possible the intrinsic properties of idealized, uniform and uncompen- 
sated (Ga,Mn)As. Our study supports the overall view of (Ga,Mn)As as a well behaved 
and understood degenerate semiconductor and band ferromagnet and, therefore, an ideal 
model system for spintronics research. We conclude in this paragraph by commenting on 
the implications of systematic studies of optimized (Ga,Mn)As materials in the context of 
the recurring alternative proposal of an intricate impurity band nature of conduction and 
magnetism of (Ga,Mn)As.'^ In the impurity band picture, the Fermi level in materials with 
~ 10^^ cm^^ Mn-acceptor densities is assumed to reside in a narrow impurity band detached 
from the valence band, i.e., the band structure keeps the form closely reminiscent of a sin- 
gle isolated Mnca impurity level. Previously, the systematic measurements of the infrared 
conductivity on the extensive set of optimized materials^^ disproved one of the founding 
elements of the impurity band picture which was the red-shift of the mid-infrared peak with 
increasing doping.f^ In the systematic measurements in Ref. HH the mid-infrared peak was 
observed to blue-shiflp2E2l experimentalists focusing on the infrared spectroscop ^^^ * ^^ *^ 
reached the consensus that the valence and impurity bands are merged in the highly doped 
ferromagnetic (Ga,Mn)As materials. The large values of the spin stiffness of the order 
meVnm^, experimentally determined in the present work, are consistent with model Hamil- 
tonian and ab initio calculation^^^'^^**^ which all consider or obtain the band structure of the 
ferromagnetic (Ga,Mn)As with merged valence and impurity bands.'^ On the other hand, 
for carriers localized in a narrow impurity band the expected spin stiffness would be small 
in a dilute moment system like (Ga,Mn)As, in which the magnetic coupling between remote 
Mn moments is mediated by the carriers.'^ By recognizing that the bands are merged, the 
distinction between a "valence" and "impurity" band picture of ferromagnetic (Ga,Mn)As 
becomes mere semantics with no fundamental physics relevance. Simultaneously, it is im- 
portant to keep in mind that the moderate acceptor binding energy of Mnca shifts the 
insulator-to-metal transition to orders of magnitude higher doping densities than in the case 
of common shallow non-magnetic acceptors J ^^ * ^^ ^ Disorder and correlation effects, therefore, 
play a comparatively more significant role in (Ga,Mn)As than in degenerate semiconductors 
with common shallow dopants and any simplified one-particle band picture of ferromagnetic 
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(Ga,Mn)As can only represent a proxy to the electronic structure of the material. 
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EXPERIMENTS PRECEEDING OUR DETERMINATION OF MICROMAGNETIC 

PARAMETERS OF (GA,MN)AS 

Magnetic anisotropy fields, Gilbert damping constant and spin-stiffness are the basic 
parameters of a ferromagnet which determine its micromagnetic properties. The anisotropy 
fields are associated with the energy required to coherently rotate magnetic moments of the 
entire ferromagnet. They can be determined in a straightforward way in magnetization or 
magneto-transport measurements from external magnetic fields required to reorient the 
magnetization of a ferromagnetic sample, or in magnetization dynamics experiments from the 
field-dependent resonant frequencies [1, 2, 3, 4]. 

Gilbert damping characterizes dissipative processes that drive the magnetization 
motion towards an equilibrium state. This phenomenon is usually investigated by the 
frequency-domain-based ferromagnetic resonance (FMR) experiment where the 
phenomenological Gilbert damping coefficient a is deduced from the resonance peak 
line width [3, 4]. The experimentally measured FMR linewidths contain not only the 
frequency-dependent linewidth due to the Gilbert damping but also the frequency- 
independent inhomogeneous linewidth broadening [3, 4]. To separate them, it is necessary to 
measure the linewidths at several microwave frequencies [3, 4]. In FMR these frequencies are 
given by the resonant-cavity frequency that significantly complicates the frequency change. 
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Therefore, the experiments are usually performed at only two different frequencies (typically, 
9 and 35 GHz [4]) that makes the corresponding separation of the individual components in 
the measured signal rather questionable. Alternatively, a can be determined from the time- 
domain based magneto-optical pump-and-probe experiment by fitting the damping of the 
measured oscillatory data by Landau-Lifshitz-Gilbert equation [5, 6]. However, to obtain the 
Gilbert damping coefficient from the measured value of a it is necessary to take into account 
a realistic magnetic anisotropy of the investigated samples (see below). Moreover, also the 
frequency dependence of a has to be measured for a separation of the intrinsic value of the 
Gilbert damping coefficient from the inhomogeneous parts of a. The absence of these two 
requirements and the un-optimized magnetic properties of the investigated samples led to a 
large scatter in the deduced values of a for Gai-xMnxAs with a different Mn content x: The 
increase of a from « 0.02 to « 0.08 for the increase of x from 3.6% to 7.5% was reported in 
Ref. 5. On the contrary, in Ref. 6 the values of a from 0.06 to 0.19 - without any apparent 
doping trend - were observed for x from 2% to 11%. 

The spin-stiffness is associated with the exchange energy of non-uniform local 
directions of the magnetization, in particular with the energy of small wave-vector spin-wave 
excitations of the ferromagnet. Considering a specific model of thermodynamic properties of 
the studied ferromagnet, the spin-stiffness can be indirectly inferred from the measured 
temperature dependence of magnetization [7], Curie temperature [7], or domain wall width 
[8]. The direct determination of the spin-stiffness from magnetization dynamics experiments 
is significantly more challenging than in the case of the magnetic anisotropy fields [9-13]. The 
low-energy non-uniform collective excitations of the system can be strongly affected by 
inhomogeneities or surface properties of the ferromagnet for which specific models have to be 
assumed in order to extract the spin-stiffness constant from the measured data. An exception 
are the Kittel spin-wave modes of a uniform thin-film ferromagnet for which the spin-stiffness 
parameter D is directly obtained from the measured resonant fields (see below). To date, spin- 
wave resonance measurements of (Ga,Mn)As have been reported on > 100 nm thick epilayers 
[9-12]. The Kittel modes with A//„ ~ were observed only in a 120 nm thick, 8% Mn doped 
(Ga,Mn)As for magnetic fields applied close to the magnetic easy-axis [11]. Measurements of 
the same sample in other field orientations showed different trends which indicated the 
presence of strong inhomogeneities and surface dependent effects [11]. A non-Kittel-like 
linear or sublinear dependence of the resonant fields on the mode index has been reported also 
in the other ferromagnetic resonance measurements of thick (Ga,Mn)As epilayers [9-12]. In 
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complementary studies of the magnetization dynamics induced and detected by magneto- 
optical pump-and-probe measurements, only two resonant frequencies were identified [13]. 
Based on the theoretical modeling, they were not ascribed to the Kittel modes but rather to 
coupled bulk-surface modes which again made the extraction of the spin-stiffness constant 
dependent on the considered model of bulk and surface properties of the studied sample [13]. 
The extracted values of the spin-stiffness from all available magnetic resonance data in 
(Ga,Mn)As materials, complemented by values inferred from magnetization and domain 
studies [7, 8], are scattered over more than an order of magnitude and show no clear trend as a 
function of Mn-doping or other material parameters of the (Ga,Mn)As ferromagnetic 
semiconductor [14]. 

In this Supplementary material we show how we are able to deduce from a single 
magneto-optical pump-and-probe experiment all these micromagnetic parameters. In 
particular, the anisotropy fields can be determined from the dependence of the precession 
frequency on the external magnetic filed combined with the probe-polarization dependence of 
the precession signal amplitude. The Gilbert damping constant can be deduced from the 
precession signal damping. Finally, the spin stiffness can be obtained from the mutual spacing 
of the precession modes which are present in the measured oscillatory magneto-optical signal. 



SAMPLES 

The time-resolved magneto-optical experiments described below were performed in a 
large set of optimized (Ga,Mn)As epilayers whose selected properties are described in detail 
in the main paper. In Fig. 1 we show results of the Hall effect measurements at 4.2 K. For this 
purpose the samples were lithographically patterned into Hall-bars of 60 ^m width. It can be 

seen in the figure that the Hall signal is affected by longitudinal magnetoresistance i?„ of the 
samples, especially at low dopings. Therefore, we extracted p from high field data and by 
fitting the measured transversal resistance Rxy by 

= B/(epd) + kiR„ + k2Rj (1) 

where d is the sample thickness and ki and k2 are fitting constants reflecting the anomalous 
Hall effect and possible imperfections in the geometry of the Hall bars. We also emphasize 
that, apart from the common experimental scatter and from the corrections due to the non-zero 
magnetoresistance and due to the anomalous Hall effect, the carrier density can in principle be 
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inferred only approximately from the slope of the Hall curve in a multi-band, spin-orbit 
coupled exchange-split system such as the (Ga,Mn)As. The error bar due to the multi-band 
nature is estimated to be -20% [15]. Due to these uncertainties we can only make semi- 
quantitative conclusions based on the measured Hall effect hole densities. 




l^nH (T) (T) 

■^0 ext ^ ' "^0 ext ^ ' 



Fig. 1. (a) Longitudinal resistances [normalized to RyJO)], and (b) transversal (Hall) resistances R„ as a 
function of normal magnetic field fioHext measured in (Ga,Mn)As epilayer with depicted Mn concentration x; 
samples temperature 4.2 K. 

For an evaluation of material parameters from an experimentally measured data (e.g., 
for an evaluation of the hole densities from the measured transversal resistances which is 
described above) it is necessary to know the (Ga,Mn)As epilayer thicknesses. However, 
accurate determination of layer thicknesses is a nontrivial task in case of thin (Ga,Mn)As 
layers. Some standard techniques (e.g.. X-ray reflectivity or optical ellipsometry) are 
inapplicable due to the weak contrast between the (Ga,Mn)As layer and the GaAs substrate, 
or unknown optical parameters. The relative accuracy of other common techniques (e.g., of 
X-ray diffraction) does not exceed 10% because of the small thickness of the measured layer. 
Therefore, we used a thickness estimation based on the following quantities: (i) the growth 
time and the growth rate of the GaAs buffer layer measured by the RHEED oscillations 
(typical accuracy of ±3%); (ii) increase in the growth rate by adding the known Mn-flux 
measured by the beam-flux monitor relatively to the Ga flux (typical accuracy of ±5% of the 
Mn vs. Ga flux ratio); (iii) reduction of thickness by the native oxidation (-1.5 nm ± 0.5 nm); 
(iv) reduction of thickness by thermal oxidation (-1.0 nm ± 0.5 nm). Relative accuracy of 
steps (i) and (ii) was verified on separate calibration growths of (Ga,Mn)As on AlAs, where 
an accurate X-ray reflectivity method to measure the (Ga,Mn)As layer thickness could be 
used. Typical thicknesses of the native and the thermal oxides in steps (iii) and (iv) were 
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determined by XPS. The resulting total accuracy of the (Ga,Mn)As layer thickness 
determination is thus 3% (relative random error) and 1 nm (systematic error). 

EXPERIMENTAL DETAILS ABOUT MAGNETO-OPTICAL EXPERIMENTS 

We investigated laser-pulse induced dynamics of magnetization by a pump-and-probe 
magneto-optical (MO) technique. A schematic diagram of the experimental set-up is shown 
in Fig. 2. The output of a femtosecond laser is divided into a strong pump pulse and a weak 
probe pulse that are focused to a same spot on the sample. Laser pulses, with the time width 
of 200 fs and the repetition rate of 82 MHz, were tuned to 1.64 eV, i.e. above the 
semiconductor band gap, in order to excite magnetization dynamics by photon absorption. 
The pump pulses were usually circularly polarized (with a helicity controlled by a wave plate) 
and the probe pulses were linearly polarized. The measured magneto-optical signals 
correspond to the probe polarization rotation induced by the pump pulses (see Fig. 2). The 
experiment was performed close to the normal incidence geometry (^, = 2° and 8° for pump 
and probe pulses, respectively) with a sample mounted in a cryostat, which was placed 
between the poles of an electromagnet. All the experimental data in this Supplementary 
material were measured at temperature of 15 K, at pump excitation intensity 30 -40 |iJ.cm'^, 
and they correspond to the helicity-independent part of the measured signal [16]. The external 
magnetic field Hext was applied in the sample plane at an angle ^^with respect to the [100] 
crystallographic direction in the sample plane (see Fig. 2). Prior to all time -resolved 
experiments, we always prepared the magnetization in a well-defined state by first applying a 
strong saturating magnetic field at an angle (pn and then reducing it to the desired magnitude 

of Hext. 



Fig. 2. Schematic diagram of the experimental set-up for a detection of the magnetization precession induced by 
an impact of the circularly polarized femtosecond laser pump pulse in (Ga,Mn)As. Rotation of the polarization 
plane of reflected linearly polarized probe pulses is measured as a function of the time delay A; between pump 
and probe pulses. The orientation of magnetization in the sample is described by the polar angle cp and azimuthal 
angle 0. The external magnetic field //„f is applied in the sample plane at an angle 
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There are several microscopic mechanisms that can lead to a precession of 
magnetization due to the impact of pump laser pulse, hi particular, very recently we reported 
on the precession of magnetization due to optical spin-transfer torque (OSTT) [16] and optical 
spin-orbital torque (OSOT) [17]. However, the most common mechanism, which is 
responsible for the oscillatory MO signals measured in the majority of (Ga,Mn)As samples at 
low excitation intensities, is the change of the sample magnetic anisotropy due to the pump- 
induced temperature increase [17] that is schematically shown in Fig. 3. Before an impact of 
the pump pulse the magnetization points to the easy axis direction [see Fig. 3(a)]. Absorption 
of the laser pulse leads to a photo-injection of electron-hole pairs. The subsequent fast 
nonradiative recombination of photo-injected electrons induces a transient increase of the 
lattice temperature (within tens of picoseconds after the impact of the pump pulse). The laser- 
induced change of the lattice temperature then leads to a change of the easy axis position [17]. 
As a result, magnetization starts to follow the easy axis shift by the precessional motion [see 
Fig. 3(b)]. Finally, dissipation of the heat and recombination of the excess holes lead to the 
return of the easy axis to the equilibrium position and the precession of magnetization is 
stopped by the Gilbert damping [see Fig. 3(c)]. The most important point from the perspective 
of the present paper is that the precession of magnetization induced by the laser pulses is 
determined by the magnetic anisotropy of the sample which makes this method an all-optical 
analog to FMR [18]. 

(a) y sample plane / (b) "7 (c) /^JH— r' "S" -y 



Fig. 3. Schematic illustration of the thermal laser pulse-induced precession of magnetization, (a) In the 
equilibrium, the magnetization points to the easy axis direction, which is located in the sample plane at azimuthal 
angle (p. (b) Impact of a pump pulse induces a transient increase of the lattice temperature that leads to a change 
of the easy axis position and, consequently, to the precession of magnetization, (c) Dissipation of the heat leads 
to the return of the easy axis to the equilibrium position. Simultaneously with this, the precession of 
magnetization is stopped by the Gilbert damping. 
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ANALYTICAL DESCRIPTION OF MAGNETIZATION DYNAMICS IN (GA,MN)AS 

The dynamics of magnetization is described by the Landau-Lifshitz-Gilbert (LLG) 
equation. We used LLG equation in spherical coordinates where the time evolution of 
magnetization magnitude and orientation, which is characterized by the polar 6 and 
azimuthal (p angles, is given by: 



dt 

d9^_ 
dt 



= , 

r 



dcp _ y 
dt "(l + a^jM^ sin^ 



V, sin 6*^ 



f 



A 

\ sin^y 



(2) 
(3) 
(4) 



where a is the Gilbert damping coefficient. The gyromagnetic ratio y = {gjUg)/h, where g = 2 
is the Lande g-factor of Mn moments, jub the Bohr magneton, and h is the reduced Planck 
constant. Functions A = rfF/rf6' and B = dF/d^ are the derivatives of the energy density 
functional F with respect to 6* and q), respectively. We expressed F in a form [2]: 



F=M 



1 . 



K.. 



-sin^ 2(p sin^ 6' + cos^ 6 -K^^,cos^ 6 ^sin^ 6'(l-sin2^)- 

v4 J 1 

■ [cos cos + sin ^ sin d^j cos{<p - <Pjj )] 



,(5) 



where Kc, and Ko^t are constants that characterize the cubic, uniaxial and out-of-plane 
magnetic anisotropy fields in (Ga,Mn)As, respectively, and Hext is the external magnetic field 

whose orientation is given by the angles 9h and q>H. For a small deviations 69 and S(p from the 
equilibrium values 6o and (po, the solution of Eqs. (3) and (4) can be written in a form 

e{t) = (9o + A^e"*'^ cos{27r ft + O J, (6) 
^t) = %+ A^e-'-^ cos(2;r/f + <D J, (7) 



where the constants A^ (A^) and 0e {0^ describe the initial amplitude and phase of 9 (<p), 
respectively. The precession frequency/and oscillation damping rate ko are given by 
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(8) 
(9) 



In our case, the investigated (Ga,Mn)As epilayers are in-plane magnets (i.e., 9 ^ n/2), the 
external magnetic field is applied in the sample plane (i.e., 6h ~ 7i/2), and the precession 
damping is relatively slow (i.e., w 0) which yields 



/ = 



^«rCOs(^-^3^)-2^:„„,+ 



kX3 + cos4(p) 



+ 2ii:„sin' 



X {fieM cos(^ -9h)+ cos 4^ - 2K^ sin(2^)) 
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(10) 



" 2h 



2H^^cos{(p-(p„)-2K^^, + ^{3 + 5cos4(p)+KXl-3sm{2(p))]. (11) 
V 2 J 



Eq. (10) express the sensitivity of the magnetization precession frequency to the magnetic 
anisotropy of the material that is a well-known effect which form the basis for the 
interpretation of FMR [3]. More interestingly, Eq. (11) shows that the precession damping ko, 
which is measured experimentally, depends not only on the Gilbert damping parametr a but 
also on the sample anisotropy and on the mutual orientation of the external magnetic field and 
the magnetization. 

We note that in previously reported magneto-optical pump-and-probe experiments [5, 
6, 19-23] the measured experimental data were modeled by LLG equation in the form 
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dM{t) 
dt 



= -r 



M{t)xH^{t) 



+ - 



a 



M{t)x 



dM(t) 
dt 



(12) 



where Heff is the effective magnetic field. However, in (Ga,Mn) As the magnetic anisotropy is 
rather complex and, therefore, modeUng of MO signals by LLG in this form does not provide 
realistic values of a because it is not possible to disentangle the effect of magnetic anisotropy 
from a [see Eq. (11)]. We believe that this is one of the reasons why the dependence of a on 
Mn concentration was so different in Ref. 5 and Ref 6. Similarly, the change of magnetic 
anisotropy of (Ga,Mn)As during the deposition of metal overlayer could be partially 
responsible for the changes of a that were reported in Ref. 23. 



EVALUATION OF MAGNETIC ANISOTROPY 

The dependence of the precession frequency on the magnetic anisotropy fields and on 
the magnitude and orientation of external magnetic field [cf. Eq. (10)] enables to evaluate the 
magnetic anisotropy from the experimentally measured precession frequencies very similarly 
as in the case of FMR [3]. In particular, for a sufficiently strong external magnetic field (p - 
(pH and the following equations can be used to fit the precession frequencies measured 

a) for Hext along the [110] crystallographic direction (i.e., q>H = 7i/4): 

/ = -2K^., + KJH^^, -2K^ -2KJ (13) 

b) for Hext along the [010] crystallographic direction (i.e., (pn - n/l): 

/ = ^ - 2^„„, + 2K^ + IH^ + 2K^ ) (14) 

c) for Hext along the [-1 10] crystallographic direction (i.e., (pH = 3n/4): 

f = ^4{H^ -2K„„, + K^+2KJH^,-2K^ + 2KJ (15) 

As an example, in Fig. 4 we show the measured dependences /(//e;^,) and their fits for two 
orientations of Hext- To increase the precision of the magnetic anisotropy determination even 
further, for all the investigated samples we supplemented this method by two additional 
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experimental techniques that provide information about the samples magnetic anisotropy - 
namely, the probe-polarization dependence of the MO precession signal amplitude and 
SQUID magnetometry. 
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Fig. 4. Dependence of the precession frequency / on external magnetic field applied along the [010] and 
[110] crystallographic directions in Gai.j.Mnj.As epilayer with x = 5.2% (points); the lines are fits by Eqs. (14) 
and (13), respectively, with A'^ = 31 mT, Ku = 27.5 mT, and A'oui = -190 mT. 

In (Ga,]V[n)As there two MO effects that are responsible for the measured rotation of 
the polarization plane zly^of the reflected linearly polarized light at normal incidence [24]. The 
first of the IVIO effects is the well-known polar Kerr effect (PKE), where A/] occurs due to the 
different index of refraction for a and a ~ circularly polarized light propagating parallel to 
the direction of magnetization M. The polarization rotation due to PKE is proportional to the 
projection of magnetization to the direction of light propagation, it is linear in magnetization 
(i.e., its sign is changed when the direction of magnetization is reversed), and it is independent 
on the orientation of the input linear polarization (see Fig. 5(b) for the angle definition) 
[24]. The second MO effect is the magnetic linear dichroism (MLD), which originates from 
different absorption (reflection) coefficient for light linearly polarized parallel and 
perpendicular to M, that occurs if the light propagates perpendicular to the direction of 
magnetization M. The polarization rotation due to IVILD is proportional to the projection of 
magnetization to the direction perpendicular to the direction of light propagation, it is 
quadratic in magnetization (i.e., its sign is not changed when the direction of magnetization is 
reversed) and it varies as sm{2fi) [24]. In Fig. 5(a) we show the MO signals measured by 
probe pulses with different orientations /] for identical pumping conditions. The measured 
dynamical IVIO signal SMO, which is a function of the time delay between pump and probe 
pulses At and the probe polarization orientation fi, can be fitted well by the phenomenological 
equation [24], 
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SMO{At,j3) = A{j3)cos[27rfAt + 0{jB)y + c{jB)e~ 



-At/ 



(16) 



where A and C are the amplitudes of the oscillatory and pulse function, respectively, / is the 
ferromagnetic moment precession frequency, is the phase factor, tg is the Gilbert damping 
time, and Tp is the pulse function decay time. All the measured data in Fig. 5(a) can be fitted 
well by Eq. (16) with a one set of parameters/, tg and Zp. The dependence obtained by 
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Fig. 5. (a) Dynamics of the MO signal measured by probe pulses with different probe polarization orientations 
in (Ga,Mn)As epilayer with x = 5.2% for jJoH^^, = mT (points); lines are fits by Eq. (16) with parameters/= 3.2 
GHz, To = 360 ps and Zp = 1050 ps. (b) Definition of the angle P that describes the orientation of the probe 
polarization plane E. 
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Fig. 6. Probe-polarization dependence of the oscillatory part A of the MO signal that was obtained by fitting the 
dynamics shown in Fig. 5(a) by Eq. (16); the values of A at time delay of 200 ps are shown (points). Lines are 
fits of A{p) by a sum of a polarization-independent signal due to PKE and a polarization-dependent signal due to 
MLD (Eq. (2) in Ref 24). The vertical arrow depicts the deduced easy axis position in the sample without the 
pump pulse, (po- 



this fitting procedure is displayed in Fig. 6. The position of the maximum in the dependence 
A{p) at 120° corresponds to the equilibrium position of the easy axis in the sample - i.e., 
the in-plane position of magnetization without the pump pulse (po [24]. The position of the 
easy axis in the sample plane is given by the relative magnitude of the cubic {K^ and uniaxial 
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(^u) anisotropy fields. Therefore, by measuring <po without external magnetic field applied, we 
are directly measuring the ratio Kc I Ku- 

The in-plane anisotropy constants can be obtained also from magnetization loops 
measured by SQUID magnetometry. For any external magnetic field the orientation of 
magnetization is determined by the minimum of the energy [cf. Eq. (5)]. If the orientation of 
magnetization as function of external magnetic field is known, the projection of the 
magnetization into the measurement axis can be easily numerically evaluated for every point 
of the magnetization loop. To obtain the anisotropy constants, we fitted the experimental data 
measured by SQUID until we obtained the best agreement between the data and the calculated 
magnetization - see Fig. 7. It is worth noting that this model does not describe the switching 
mechanism (governed by the domain wall physics which is not treated in our single domain 
description), so the parts of hysteresis loops containing the switching were not used in the 
analysis. Moreover, in the case of uniaxial systems (Ku>Kc) an analytical expression for the 
magnetization measured along the hard axis can be utilized to analyze the data [1]. 
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Fig. 7. Evaluation of the magnetic anisotropy from the SQUID magnetometry. The SQUID measurement along 
[110] crystallographic direction in (Ga,Mn)As epilayer with x - 5.2% (points) is compared with the calculated 
magnetization projection for A'^. = 31 mT and A'„ = 27.5 mT (line). 

To sum up, by a simultaneous fitting of the measured dependence of the precession 
frequency on an external magnetic field (Fig. 4), of the MO signal precession amplitude on a 
probe-polarization (Fig. 6), and of the data measured by SQUID magnetometry (Fig. 7) we 
evaluated very precisely the magnetic anisotropy for all the investigated samples. The 
example of the obtained in-plane angular dependence of the free energy in (Ga,Mn)As 
epilayer with Mn concentration x = 5.2% is shown in Fig. 8. 
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DETERMINATION OF GILBERT DAMPING COEFFICIENT 

For numerical modeling of the measured MO data, we first computed from the LLG 
equation (Eqs. (3) and (4) with the measured magnetic anisotropy fileds) the time-dependent 
deviations of the spherical angles [S6(t) and S(p(t)] from the corresponding equilibrium values 
{do, (Po)- Then we calculated how such changes of 9 and (p modify the static magneto-optical 
response of the sample MO^"", which is the signal that we detect experimentally [24] : 

^0(a?,^) = -^(aOp^^^ + ^(aOp'''^2cos2(^, -^)+^^^^P^^^2sin2(^?)o -^). (17) 

The first two terms in Eq. (17) are connected with the out-of-plane and in-plane movement of 
magnetization, and the last term describes a change of the static magneto-optical response of 
the sample due to the laser-pulses induced demagnetization [24]. p^^^ and p^^'^ are MO 
coefficients that describe the MO response of the sample which we measured independently 
in a static MO experiment for all the samples - see Fig. 9 for MO spectra measured in sample 
with x = 5.2%. 
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Fig. 9. Spectral dependence of static PKE and MLD in (Ga,Mn)As epilayer with x = 5.2%, the arrow indicate the 
spectral position of the laser pulses used in the time-resolved experiment shown in Fig. 5 and Fig. 10; note that 
the data for MLD are multiplied by 5 for clarity. 

The examples of the fitting of the dynamical MO optical data are shown in Fig. 10. 
The measured data can be fitted well by LLG for time delays longer than -150 ps, which is a 
time that it takes to establish the quasi-equilibrium conditions in the sample. We stress that 
the only fitting parameters in our modeling are the Gilbert damping coefficient a, the initial 
deviation of the spherical angles from the corresponding equilibrium values, and the 
parameters describing the in-plane movement of the easy axis and the demagnetization signal, 
which are apparent as the non-oscillatory signal in the measured dynamics [24]. The obtained 
dependence of a on Hext is shown in Fig. 1 1(a) for two different orientations of Hext- For Hexi 
applied along the [010] direction, a decreases monotonously with Hext- On the contrary, for 
Hext applied along [1 10] direction, a is a non-monotonous function of Hext reaching a similar 
values of a for mT and 100 mT. However, this non-monotonous dependence is a 
consequence of the field-induced frequency decrease (see Fig. 4) when the magnetic field is 
applied along the magnetically hard [110] direction (see Fig. 8). When a is plotted as a 
function of the precession frequency (rather than the external field) we do not observe any 
significant difference between the different crystallographic directions - see Fig. 11(b). A 
field dependent damping parameter was reported in various magnetic materials and a variety 
of underlying mechanisms responsible for it were suggested as an explanation [25-29]. We 
note that the damping parameter a extracted from the fits should be regarded as a 
phenomenological parameter that accounts for combined effects of a (frequency independent) 
intrinsic Gilbert damping, an inhomogeneous broadening, a two magnon scattering, and 
various propagation spin wave processes resulting from the nonuniform spatial profile of the 
excited precession. We also note that the rate of decrease of a with/ is sample dependent and, 
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therefore, we show in the main paper the doping dependence of the frequency-independent 
part of a. 




Time delay (ps) Time delay (ps) 

Fig. 10. Dynamics of the MO signal measured for external magnetic field (a) jJfflea - 10 mT and (b) iMflext = 
400 mT applied along the [010] crystallographic direction in (Ga,Mn)As epilayer x = 5.2% (points); lines are fits 
by LLG. 




Fig. 11. (a) Dependence of the Gilbert damping coefficient a on external magnetic field H„, applied along the 
[010] and [1 10] crystallographic directions in (Ga,Mn)As epilayer with x = 5.2%. (b) Same data as in (a) but as a 
function of the precession frequency/. 



DETERMINATION OF SPIN STIFFNESS 

As we show in Fig. 5 of the main paper, we observed more than one precession mode 
in (Ga,]VIn)As epilayers with a sufficient thickness. These precession modes are the spin 
wave resonances (SWRs) - i.e., spin waves (or magnons) that are selectively amplified by 
fulfilling the boundary conditions of the thin magnetic film [18, 30]. Up to now, SWRs in 
(Ga,]V[n)As were investigated mainly in a frequency-domain where they are apparent as 
multiple absorption peaks in the FIVIR spectra [3, 10 - 12]. The existence of multiple 
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resonances in FMR reveal that there exist several external magnetic fields at which the 
Larmor precession frequency in the sample coincides with the microwave frequency. The 
resonant field for the n-th mode (Hn) is obtained by solving the LLG equation with a term 
corresponding to exchange interactions in the material and by considering the appropriate 
boundary condition [11]. In homogeneous thin films with a thickness L, only the 
perpendicular standing waves with a wave vector k fulfilling the resonant condition 
kL=n7rsiTe amplified; the mode with n = denotes the uniform magnetization precession 
with zero k vector. In principle, there exist two symmetric boundary conditions which are 
schematically illustrated in Fig. 12. The position of n-th SWR mode in the FMR spectrum H„ 
is given by the Kittel relation [11] and the following equation applies 

^n^H,-H,^=n'^'^, (18) 

where n is an integer, D is the exchange spin stiffness constant, jUb is the Bohr magneton, g is 
the g-factor, and L is the sample thickness 



(a)n=0n=1n = 2 ^^^^^q^^^^^2 




Fig. 12. Spin wave resonances in homogeneous thin magnetic films with a thickness L that have a node (a) or 
maximum (b) at the surface; n is the mode number. 

In FMR only the modes with odd n are observed [11] and the corresponding resonant 
fields are smaller than that of the uniform magnetization precession (i.e., AH„ > 0). In the 
magneto-optical pump-and-probe experiment, the external magnetic field is kept constant 
during the measurement of any dynamical MO trace. Consequently, the SWRs are apparent as 
additional frequencies that are larger than that of the uniform magnetization precession. 
Ultrafast optical pulses also excite all resonant modes without any k selectivity [18, 30]. 
Consequently, for a homogeneous magnetic film with a given thickness, a higher number of 
SWRs is detectable in the MO dynamical traces than in the FMR spectra. This is particularly 
important for (Ga,Mn)As that is magnetically homogeneous only when prepared in a form of 
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rather thin films and, therefore, where only a limited number of SWRs is present within a 
detectable range of the precession frequencies. For an external magnetic field Hgxt appUed in 
the sample plane, the angular frequency of the n-th SWR mode/„ is given by [3, 31] 



f _ 

■'" h 



i(//^cos(^-^J-2/..,. ^-(^"^-"-^^) .2/.„sm-[^-^].Mf; 
X (^^ cos(^ - g)ij)+2K^ cos 4^ - 2K^ sm{2^)+AH„ ) 

(19) 
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which enables to convert the experimentally measured frequency spacing of individual modes 
to the field differences AH„ from which the magnitude of the spin stiffness D can be evaluated 
using Eq. (18) (see Fig. 5 in the main paper). 

As we illustrate in the following chapter, the magnetic homogeneity of the 
investigated epilayer is absolutely essential for a correct determination of D from the 
measured SWR spacing. Therefore, the experimental results obtained in samples that had 
been prepared by etching the original 48 rmi thick (Ga,Mn)As epilayer down to the thickness 
39, 29 and 15 nm are of fundamental importance. In Fig. 13 we show the corresponding FFT 
spectra of the measured oscillatory MO signal. Clearly, the frequency fo of the lowest SWR 
does not depend on the film thickness. This confirms that the lowest observed SWR really 
corresponds to the uniform precession of magnetization and, moreover, it proves that this film 
is magnetically homogeneous. Also the spacing AHj shows the expected [see Eq. (2)] linear 
dependence on and 1/L^ (see Fig. 5 in the main paper) that enables a reliable determination 
of the value of D. In Fig. 14 we compare the experimental data for 48 nm and 15 nm thick 
epilayers from which the FFT spectra depicted in Fig. 13 were computed. Clearly, the etching 
of the sample from 48 nm to 15 rmi not only suppressed the higher SWRs, which is apparent 
from the purely sinusoidal shape of the data for the 15 nm fihn, but it also increased the 
precession damping, which is probably a consequence of a slight variation of the etched fihn 
thickness within the laser spot size of 25 (im. These data illustrate that the magneto-optical 
pump-and-probe experiment is a very sensitive diagnostic tool not only of the magnetic but 
also of the structural quality of thin magnetic fihns. 
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Fig. 13. Fourier spectra of oscillatory MO signals measured for iMMext = 20 mT applied along the [010] 
crystallographic direction in samples prepared by etching from 48 nm thick (Ga,Mn)As epilayer with x = 7% 
(points), the curves are labeled by the film thicknesses, normalized and vertically shifted for clarity; the lines are 
fits by a sum of Lorentzian peaks. 
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Fig. 14. Comparison of oscillatory parts of MO signals measured in the original 48 nm thick epilayer and in the 
epilayer that was etched down to 15 nm; the curves are normalized and vertically shifted for clarity. 
Experimental conditions are described in Fig. 13. 



DEMONSTRATION OF INAPPLICABILITY OF SPIN STIFFNESS 
MEASUREMENT IN THICK (GA,MN)AS EPILAYERS 

Finally, we illustrate the significance of the film magnetic homogeneity for a correct 
evaluation of the spin stiffness. For this purpose we selected a 500 nm thick (Ga,Mn)As 
epilayer with 7% Mn (i.e, a sample with the same nominal Mn doping as the one used in 
experiments depicted in Fig. 13 and Fig. 14). In Fig. 15 we show the temperature dependent 
magnetization projections to several crystallographic directions measured in the as-grown and 
annealed samples. In the as-grown sample, the temperature dependence of magnetization 
projections is strongly non-monotonous [see Fig. 15(a)]. Moreover, the Curie temperature Tc 




18 



is only ~ 60 K that is very low for a material with 7% Mn. This is a consequence of a high 
concentration of unintentional interstitial Mn impurities in the sample that compensate both 
the local moment and the holes produced by substitutional Mn atoms [32]. The amount of 
interstitial Mn impurities in the sample can be reduced by a thermal annealing [33]. However, 
even very long annealing times are not sufficient for obtaining a high quality sample from the 
thick epilayer due to the formation of the surface oxide that controls the outdiffusion of 
interstitial Mn impurities [33]. Therefore, the 40 h long annealing at 200 °C led to an increase 
of Tc but only to 90 K, which is still substantially lower than Tc 150 K observed in thin 
samples with the same nominal concentration of Mn. Simultaneously, the temperature 
dependence of magnetization does not show the expected sharply vanishing magnetization at 
Tc (cf. Fig. 2 in the main paper for the data in optimized epilayers). 
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Fig. 15. Temperature dependence of the magnetization projections to different crystallographic directions 
measured by SQUID in 500 nm thick (Ga,Mn)As epilayer with x - 1%. (a) As-grown sample, (b) Sample 
annealed for 40 hours at 200°C. 
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Fig. 16. Temperature dependence of the resistivity p (a) and its temperature derivative dp/dT (b) measured in 
500 nm thick (Ga,Mn)As epilayer with x = 7% annealed for 40 hours at 200°C. 



In Fig. 16 we show the temperature dependence of the resistivity and its temperature 
derivative measured in the annealed sample. Clearly, there is no sharp Curie point singularity 
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in the temperature derivative of the resistivity which is the fingerprint of a high magnetic 
quality of (Ga,Mn)As epilayer (cf. Fig. 1(a) in the main paper) 

In Fig. 17 we show the time-resolved magneto-optical signals measured in this 500 nm 
thick epilayer. In the as-grown sample two precession modes can be identified. In the 
annealed sample the improved magnetic quality leads to a strong suppression of the 
magnetization precession damping with respect to that observed in the as-grown sample. For 
example, the data shown in Fig. 17(a) and (b) for the lowest modes correspond to damping 
times of 210 ps and 460 ps for the as-grown and annealed sample, respectively. In addition, 
the annealing led to a considerable increase of the number of observed SWR modes in the 
measured TRMO signal. However, their identification is a rather complicated task. 




-60-U ' 1 ' 1 U I 

500 1000 1000 2000 3000 

Time delay (ps) Time delay (ps) 

Fig. 17. Time-resolved magneto-optical signals (points) measured in as-grown (a) and annealed (b) 500 nm thick 
(Ga,Mn)As epilayer with x = 7%; note the different jc-scales in (a) and (b). The lines in the upper parts of the 
figures are a sum of damped harmonic functions and the corresponding precession modes are plotted in the lower 
parts of the figures. External magnetic field of 10 mT was applied along the [010] crystallographic direction. 

In Fig. 18 we show the FFT spectrum of the oscillatory MO signals measured in the 
annealed sample for external magnetic fields of 10 mT and 20 mT. Even though the magnetic 
field change was rather small, the FFT spectra were changed dramatically. In particular, at 10 
mT there are 3 peaks with comparable intensities (and 5 peaks in total) while at 20 mT there 
is only 1 strong peak (and 4 peaks in total). In Fig. 19 we show the dependence of the 
frequency of SWR modes on the external magnetic field - at least for the first sight, it is not 
apparent how to assign the observed modes to mode numbers defined by Eq. (18), which is 
the basic requirement for an evaluation of the spin stiffness from the measured data. 



20 




Fig. 18. Fourier spectrum of the oscillatory part of the MO signal measured in the annealed sample for external 
magnetic fields of 10 mT (a) and 20 mT (b) applied along the [010] crystallographic direction (points); the red 
line is a fit by a sum of Lorentzian peaks (green lines) and the arrows indicate positions of the peak frequencies. 
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Fig. 19. Dependence of the precession frequency /„ on measured in the annealed sample for external 
magnetic field applied along the [010] crystallographic direction. 



In Fig. 20 (a) we show a plausible assignment of the measured frequencies to four 
SWRs described in the previous chapter and one non-propagating surface mode [11]. We note 
that the identification of the lowest mode for fields below 15 mT as the surface mode is based 
on the analysis reported in Ref. 1 1 - in particular, due to the observations that this mode is 
apparent only at certain external magnetic fields and that it has a smaller amplitude that the 
one assigned to the homogeneous precession [see Fig. 18(a)]. Following the analysis reported 
in the previous chapter, we can now proceed to the evaluation of the spin stiffness. In Fig. 
20 (b) the deduced values of zlf/„ are plotted as a function of n . The observed mode spacing 
deviates significantly from that expected for SWRs in a magnetically homogeneous film [see 
Eq. (18)] which is another fingerprint of the magnetic inhomogeneity in this 500 nm thick 
epilayer [9-12]. Consequently, despite a large number of SWRs detected in this sample, they 
cannot be used for a direct determination of the spin stiffness. 
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Fig. 20. (a) Dependence of the precession frequency /„ on //„, measured in the annealed sample for external 
magnetic field applied along the [010] crystallographic direction re-plotted from Fig. 19 with the depicted 
assignment of precession frequencies to the individual SWRs and to the surface mode, SM (points). Lines are fits 
by Eq. (19). (b) Dependence of mode spacing A//„ on square of the mode number n (points), line is the 
theoretical dependence AH„ ~ n. 
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